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We calculate the plasma excitations in a bundle as well as a two-dimensional (2D) periodic 
array of aligned parallel multishell nanotubes on a substrate. The carbon nanotubes are oriented 
perpendicular to the substrate. The model we use for the system is an electron gas confined to 
the surface of an infinitely long cylinder embedded in a background dielectric medium. Electron 
tunneling between individual tubules is neglected. We include the Coulomb interaction between 
electrons on the same tubule and on different tubules for the same nanotube and neighboring 
nanotubes. We present a self-consistent field theory for the dispersion equation for intrasubband and 
intersubband plasmon excitations. For both the bundle and 2D array of aligned parallel nanotubes, 
the dispersion relation of the collective modes is determined by a three-dimensional wave vector with 
components in the direction of the nanotube axes and in the transverse directions. The dispersion 
equation is solved numerically for a singlewall nanotube 2D array as well as a bundle, and the 
plasmon excitation energies are obtained as a function of wave vector. The intertube Coulomb 
interaction couples plasmons with different angular momenta m in individual nanotubes, lifting the 
±m degeneracy of the single-nanotube modes. This effect is analyzed numerically as a function 
of the separation between the tubules. We show that the translational symmetry of the lattice is 
maintained in the plasmon spectrum for the periodic array, and the plasmon energies have a periodic 
dependence on the transverse wave vector qx. For the bundle, the Coulomb interaction between 
nanotubes gives rise to optical plasmon excitations. 

PACS numbers: 73.20.Mf, 73.61.Wp, 71.45. Gm, 61.46.+W 

I. INTRODUCTION 

Carbon nanotubes (CNTs) are an allotrope of carbon which have been prepared in various configurations. For 
example, a one-atom thick sheet of graphene rolled up into a seamless cylinder with diameter on the order of a 
nanometer constitutes a singlewall carbon nanotube (SWNT). The ratio of the length of the resulting nanostructure 
to the diameter exceeds 10^. The novel properties of such cylindrical carbon tubules have made them potentially 
useful in many applications in nanoscience. These include photonics, electronics (electrical circuits) and other areas 
of materials science (to strengthen polymer materials) Their extraordinary strength, unique electrical properties 
and their ability to conduct heat efficiently have generated considerable interest among both experimentalists and 
theoreticians.'^'^ Nanotubes may also be produced by synthesis methods, for instance multi-wall nanotubes (MWNTs) 
which consist of multiple layers of graphene rolled in on themselves to form a cylindrical tube. There are two models 
which can be used to describe the structures of multi-wall nanotubes. In the first model, sheets of graphene are 
arranged in concentric cylinders, e.g. a SWNT nanotube within a larger SWNT nanotube. In the Parchment model, 
a single sheet of graphene is rolled in around itself. The interlayer distance in multi-wall nanotubes is close to the 
distance between graphene layers in graphite, approximately 3.3 A. Double-wall carbon nanotubes (DWNT) are 
particularly interesting because they com bin e similar morphology and properties compared to SWNT, but improving 
significantly their resistance to chemicals.'-^ It is abundantly clear that carbon nanotubes are unique nanostructures 
with remarkable mechanical and electronic properties.'^l^ Interest has focused on them as prototypes for a one- 
dimensional quantum wire as well as how these one-dimensional properties would be modified when the nanotubes 
are combined in a linear array on a two-dimensional (2D) plane or arranged on a three-dimensional lattice. In Fig. 
[T] we show schematically an array of nanotubes. We show that the system of nanotubes would lead to a coupling 
between the plasmon modes. 

Over the years, there have been a number of experimental measurementd^iHlD g^j^^j^ theoretical calculationd^SHSSI 
dealing with the effects due to coupling between nanotubes. In the work of Kodak, et al.!^ employing electron energy 
loss spectroscopy (EELS), it was demonstrated that due to the strong intra-tube and inter-tube Coulomb coupling, 
there is an in-plane mode and an out-of-plane mode for an array of nanotubes arranged on a two-dimensional surface. 
This is in agreement with the work of Gumbs and AizirP who derived the dispersion equation for nanotubes whose 
axes are aligned on a planar surface (see also the paper by Shyu and LirP^. Stephan, et al.'^^ reported on a detailed 
experimental investigation the plasmon excitations of multi-wall carbon nanotubes using EELS. The results in ReiW^ 
compared the plasmon modes for MWNTs with those measured for SWNTs. In interpreting their datcP^, a continuum 
dielectric model was adopted in which the bulk response function for the material making up the nano-particles was 
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used. The EELS experiments yielded high energy excitations with energies in the 5-15 eV range. Lower energy 
plasmon excitations 0.5 eV) were measured by Murakami, et al.'^Slin ahgned SWNTs on a planar surface forming 
a thin film by means of optical absorption techniques, thereby demonstrating the anisotropy of the mode^. Films of 
thickness 1000 A consisting of SWNs with mean diameter 4 nm have been probed using EELS and the plasmon mode 
dispersion was obtained^^. This means that there are about twenty-five SWNs stacked on top of each other. Within 
the limits of the resolution for the momentum and energy transfer for the EELS experiments reported in Ref.'^. 
(See also the paper by Thess, et al.^^.) The observed plasmons are in the energy range of 5-7 eV and wave number 
~ 0.15 and are believed to be excitations of the 7r-electrons which are formed by the 2pz orbitalJ^^. For metallic 
and narrow-gap semiconducting SWNTs, the plasmon excitation energies are of low frequency with energies in the 
range < 1 eV^. These collective excitations are due principally by exciting the charge carriers in the low-energy 
bands near the Fermi level. With higher spectroscopic resolution, these low-frequency plasmon excitations should be 
observed, as we demonstrate in this paper with the use of a simple model which is valid for this range of frequencies. 
The results for the low-frequency plasmon excitations in the tight-binding approximation agree with the electron gas 
model for SWNTs and DWNTs (see Refs.^^^ ancP^ as well as references therein) and this is the reason we use the 
present model for an array. 

The electron gas model for metallic carbon nanotube bundle was presented by Lin, et alW^ following the reported 
experiments of Thess, et al.^^. Their formalism for calculating the dispersion equation for this three-dimensional 
electron system was based on the RPA. However, since the plasmon excitations cannot be categorized as intrasubband 
and intersubband plasmons as for a single tubule, it is crucial to include the coupling between subbands with different 
angular momentum quantum numbers. It is the purpose of this paper to calculate the dispersion relation of the low- 
frequency plasmons for dopecP^ nanotubes in which the charge carriers are introduced onto the graphene tubules by 
means of intercalation, as can be done in carbon fibers or CgcP^- Also, we investigate the way in which this dispersion 
relation is affected by their geometrical arrangement, such as their separation and angular configuration. This effect 
may be analyzed through the behavior of the collective plasma excitations. We do so by calculating numerically the 
plasma dispersion as a function of the separation between the tubules. We demonstrate that for a bulk lattice the 
translational symmetry of the lattice is maintained in the plasmon spectrum, and the plasmon energies have a periodic 
dependence on the transverse wave vector components. 

We first calculate the plasmon dispersion equation for an isosceles right triangle of nanotubes as a simple example 
of a bundle of nanotubes on a substrate with their axes aligned parallel to the z axis.'^^'^^H xhe weak van der Waals 
interaction which holds them together is neglected in our calculations. The carbon nanotubes are oriented perpendic- 
ular to the substrate. With one of the axes of the tubules on the z axis, the other tubules have their axes at x = a^; 
on the X axis and y = on the y axis. In the limit Oa, — >■ oo or Oj, — > oo, our dispersion equation reduces to that 
for a pair of multishell nanotubes, as obtained previously.!^ We examine how the Coulomb interaction effects on the 
plasmon spectrum depend on the axial separation. In order to consider the role played by the periodicity of a lattice 
on the collective modes, we consider an array of nanotubes, with their axes parallel to the z-axis and equally spaced 
by distance in the x-direction and ay in the y-direction. Each nanotube may consist of M > 1 co-axial cylindrical 
tubules. We assume that there are no electrons tunneling between the tubules in each nanotube and between the 
nanotubes. The plasmons are determined by the angular momentum quantum number m corresponding to transitions 
within a subband {m = 0) or between different subbands (m 7^ 0) as well as the wave vector along the axis of the 
nanotube and the transverse wave vector components and qy. 

Even in the absence of tunneling, for finite and ay, the plasmon spectrum for one nanotube is modified by the 
Coulomb interaction between the nanotubes.l^ Furthermore, the inter-tubule Coulomb interaction causes the angular 
momentum not to be conserved in a nanotube array or bundle, and modes with different m are coupled to one another. 
In particular, the degeneracy of the modes with angular momentum quantum number m and —m is lifted by the 
Coulomb interaction. 

The numerical calculations we present were carried out to analyze the way in which mode-coupling modifies the 
plasmon spectrum. For the bundle consisting of three nanotubes whose axes are at the vertices of a triangle, we 
calculate the dispersion relation as a function of the transverse and longitudinal components of the wave vector. 
For the array of nanotubes on a 2D lattice, wc show that the symmetry of the lattice is maintained in the plasmon 
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spectrum, and the plasmon excitations depend on the wave vector qj. and qy with period — and — in the x and 

ax ay 

y-directions, respectively. Numerical results are presented for the array of nanotubes on a 2D lattice with Af = 1. We 
obtain data for the plasmon excitation energies as functions of qz, q^ and the lattice period a^- Different plasmon 
modes associated with the intrasubband and intersubband electron transitions in the tubules are considered and 
compared with the results for a single nanotube and a linear array on a 2D planar sheet. 

The outline of the rest of this paper is as follows. In Sec. |llj we derive the dispersion relation for plasmons in 
a bundle of nanotubes In Sec. |III[ we present and discuss the plasma dispersion equation for a periodic array of 
nanotubes consisting of an arbitrary number of concentric tubules. In Sec. |IV[ we present numerical results for the 
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plasmon excitation energies and discuss their properties as functions of their wave vectors and the separation of the 
nanotubes. We give a summary in Sec. [Vj 

II. THEORETICAL FORMALISM FOR A BUNDLE OF TUBULES 

We first consider a simple model of bulk quantities of nanotubes. For this we have a coupled triad of infinitesinially 
thin nanotubules, with their axes aligned in the z-direction. The nanotubes are predicted to be semiconducting or 
metallic depending on the chirality of the tubules. The axis of one of the tubules of radius Ri is located at the origin 
at a; = and the remaining two tubules each of radius R2 and i?3 are located at a distance of a; = a^; and y = Uy from 
the origin on the x and y-axis, respectively. We shall impose the condition that Gx > Ri + R2 and ay > Ri + R3. 
There is no tunneling between the tubules so that the eigenfunctions of an electron on the j*'' tubule (j = 1,2) on 
the 2:-axis or y-axis , with axial wave vector and angular momentum quantum number are given by 

Wd)y = (P- , (1) 



with v = {fcz, 1} and 



"^iAp) = ^^'^0 (P) , = 5{P - Rj) ■ (2) 



The energy eigenvalues are 

The plasmon dispersion equation can be obtained by solving the density matrix equation ih— = \H,p] where H = 

at 

Ho - e0oxt and p ^ po + Sp with {vj\Ho\v'j') = ^v.jS^^'Sjj^ and {vj\po\v'j') = fo{e^j)S^^'Sjj' , where fo{eu,j) is the 
Fermi-Dirac distribution function. In the lowest order of perturbation theory, we obtain 

{v3\5p,Wf) = 2e ^f''^''^-^f''^'\ ,.Mtot{r,u^)W'f} . (4) 
huj - e^j + ejyj- 

Here, (/)tot(r, = (/)ind(r,a;) + 0cxt(r, w) is the sum of the external and induced potentials with 0ind(r, a solution 
of Poisson's equation 

4:7r6 

V 0i„d(r,w) = (5nind(r,w) . (5) 

Now, we obtain the following matrix elements for the tubules of radii R2 and i?3, i.e., 



(6) 



where q± = \ (1^ + qi- This yields 



2e 



J, 3 v,v' 



^^(z.',/|e-''i->,j)</>i„d(q,u;)( 



q,q' 



V ,J ) 



(7) 



and 
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X] '/'md('?i:,9y,gz;t^)Jm('7li?l) 
2e /o(efc^,i;2) — /o(efej-g^,Z-m;2) 



9=, 
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— ■ Fi{m,q^,^2;R2) 

X ^ 0i„d(qi:,gy,gz;w)e""'-''="i^i(rn,ql,'!?2;i?2) 

X ^ 0i„d(g;,g;,g.;w)e<'^''i^2(m,(7l,z93;i?3) , (8) 

where ^2 is the angle between the a;-axis and a hne drawn from the center of the cyhnder at the origin to a point 
on circumference of the cyhnder with its center at a; = Oa^. We also define i?3 as the angle between the y-axis and a 
line drawn from the center of the cylinder at the origin to a point on the circumference of the cylinder located on the 
j/-axis at y = ay. In addition, we have introduced the following notation: 

Fi{m,q^,d2;R2) = e^^^V^Cq^i^s) 

F2im,q^,^3;R3) = e™''V™(Qj.i?3) . (9) 
The potential 0ind(q, w) in Eq. Q can be rewritten in terms of Sniadi<l,^)- This gives 

Suindiq.,^) = - —Xl,m{Qz,i^)J7n{q±Rl)Ui,jn{qz) [ ^^^1^ 

£s \ q± 

- X2,m{lz,uj)e "^^"^ Fi{m,q±,'d2; R2)U2,m{qz) 

- —X3,m{(lz,^y-'''''''F2{m,q^,^3;R3)U3,M, (10) 



with 



, . 1 V- Snind{qx,qy,qz) . / p n <lx - i qy 

Ul,m{qz) = -j—r 1^ „2 ,„2 , „2 Jra[q^Rl) I — — 
^xi^y ^ „ qx + % + 1z \ (11- 



LxLy 



^nind{qx, 


qy,qz) 


Qx + qI 




Sriiridiqx, 


qy,qz) 


ql+ql 


+ ql 


5niad{qx, 


qy,qz) 



TT I \ ^ dnind[qx,qy,qz) iq^a^ / Q D \ 

U2,m{qz) - -j—r 2^ „2 ,„2 ,„2 ^ ^ ^^1 K 9^ , 1?2 ; i?2 ) 



C/3,™fe) - tV E '^Ti'^J e-^^-^f2(m, g^, ^a; R,) (11) 



-^^^2/ o o 9a + 9y + Iz 



and 



After substituting, we obtain 
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2e' 



{qz,^)Ui^m{qz)-j—j- XI ~J 



1 Jm'iq'MJMRi) (q'x + n 



ixiy^f^, q'x^ + q'y^ + ql 



/ \ m —m 
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2e' 



{qz,^)U2,r,i'{qz)j^ 2^ 2 4_ „2 



qx+qt 



q-x + iq^ 
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2e' 



-X3, 



1 

{qz,^^)U3,,n'{qz)-r-r X 



' " q'.^ + ql \ q'i_ 



y 



(13) 



1 . , 

U2,milz) = Xl,7n{qz,i^)Ul.m{qz)-i—r- 



,^Fi(m,9l,t92;i?2)Jm'(gli?i) 



qx + 



9± 



2e- 



2e2 



1 

-Xz..m'{qz,uj)Us.^m'{qz)-,—r- V 



(m' , gl , , i?3 ) J^i (™, gl , 1^2 , i?2 ) 



91' + 9? 



(14) 



1 . , 

t^3,m(g^) = Xl,m'{qz,^)Ui^rn'{qz)-j—j- 



,^ F2{m,q'^,^3;R3)JnAq'±Ri) 
q'± + ql 



qx + ^iy 



2e2 

El 

2e2 
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'{qz,oj)U2,„i'iqz)-j-^ 



FAm',q'^,d2,R2)F2{m,q'^,th,R3) 



LxLy 



9l' + 9.2 



X3w(9.,C.)f/3,„.(,.)-V E ^2*("^^-^^'^3,i?3)f2(m,,l,^3,i^3) ^-..,^..^.,,.. 
J^xJ^y I I 



9l' + 9.2 



(15) 



We note that the set of equations in Eqs. ( 13 1 - ( 15 1 form a a simuhaneous system in the variables Ui^m{qz)i U2.m{qz) 



and Us^miqz), with m — 0, ±1, ±2, • • •. These equations may be simpUfied if we now introduce the fohowing quantities 
below, i.e., if we set 



A 



B„ 



2e2 1 ^ J,n'{q'±Rl)Jrn{q'±Rl) f qx+^q'y'^"' 

-Xl.m'[qz,i^)-j—j— > - 
J-ixi-iy I , 



2e2 



9l' + qi 
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2e2 



-X2 



. ,x 1 ^,,^,^ fl(m^gl,^92,il;2)^^(gli?l) ( q'x + n 



V 



-X3 



. , 1 _.,'a,^2(™',gl,l93,il'3)Jw(gli?l) ^9^+^^ 

,m'[qz,^) y y e y - 



LxLy 



Ql' + q'z 



9± 
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H„ 



2e2 
2e2 



X2,m'(9z>^ 
X3,m'(9z>^^ 



1 Fi*(m',gl,7?2,i?2)i^i(m,(7^,z?2,i?: 



LxLy 



E 



+ 9? 



«'.2 _! „2 



^^^a ^r^, + gj 



LxLy 



i^2 (™, g^ , , i?3 ) ^m' (gl ^1 ) / + ^9: 



gl' + g.^ 



F^{m,q'^,d2,R2)F2{m,q'^,^3,R3) 



LrpL. 



■zi.?; 



gl'+g.^ 



gl 



1 -F2*(m', gl, t^3, R3)F2{m, g^, 1^3, -R3) -ig„a„ »?'a„ 



(16) 



We note that Fmm and Hmm are the only matrix elements containing a product of the two factors exjp{iq±ax cos'd) 
and exp(— igj^Oj^ sini?). The Jacobi- Auger formula along with a standard formula for exp(i2:cosi9) and exp(^isin??) 
produces the above obtained results. 



Back-substituting the results for the matrix elements in Eq. (16 1 into Eqs. (13), (14) and (15), respectively, yields 
the following compact system of simultaneous linear equations 

m' m' m' 

U2.m + Dfnm'Ul,m' + -E'mm'C^2,m' + Fmm'U^^m' = 

7n' 771 ' 777' 

C^3,m + Gmm'Ul^m' + ^^712777' C^2,7ra' + Imm'U3,m' — . (17) 
771' 777' 777' 

The mode coupling arising from the Coulomb interaction is clearly seen through the finite values of the coefficients 
when to' ^ m in Eq. ( |17[ ). The dimension of the associated coefficient matrix above is dependent on the type of 
transitions being investigated. Therefore, to obtain a matrix of modest dimensions, let us consider only intrasubband 
transitions (i.e., m = 0). For such transitions, we have 



/C/i,o(g^)\ /I + ^0,0 Bofi Co.o \ /C/i,o(gz 

^(W,g.,g7;,g.) U2.o{qz) - ^0,0 1+^0.0 ^^0,0 [/2,o(g.) ) = ( ) . (18) 

\(73.o(g.)/ V Go,o Ho,o l + h.J \U3fl{q, 




For non trivial solutions of Eq. ( 18 ) to exist, one requires that the determinant of the coefficient matrix e (g^,, , g^) 
vanish. That is Det ti^qx, qy, qz) = 0. This is the dispersion formula for the plasmons and particle-hole mode 
intrasubband transitions only. Let us further investigate the intrasubband (to' = to = 0) plasma dispersion equation. 
In general, an explicit determination of the elements of the matrix yields for to = m' leads to a simplification of the 



results in Eq. ( 16 ) and we have 



^777777 = l7n{qzRl)KrniqzRl)Xl.rn{qz,^) 



g2 poo 

B,nm = X2.7r7(g2,w) / dq±_ q±_ 

T^Ss Jo 



Jo{q±ax)Jm{q±Rl)Jrn{q±R2) 



il + q^z 



^ , . f°° , Joiq±a,y)J,n{q_i_Ri)J,niq±R3 

6,77777 = X3, 777(g2,^^) / C(g± g± 



TTEs 'Jo Q1+ Tz 

, \ /"°° J Joiq±ax)Jmiq±Ri)Jrn{q±R2) 



Dm-m = Xi,777(g2,w) / dq± g_L- 

TTEs Jo 



qi + Qz 
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Emm = Ir7i{qzR2)K,n{qzR2)X2,m{lz^^) 



e 



2 



oo 



i''Jk{q±ax)Jk'iq±ay)Jm{q±R2)Jmiq±R3) 



^ , s f°° , Jo{q±ay)Jm{q±Ri)J,n{q±R3) 

Ti'es Jo qi + qi 

(? ^ '^^{-^f' Jk{qi^a^)Jk'{q^av)Jrn{q^R2)Jm{qi.R-i) 

-fmm = Im{qzR3)K„i{qzR5)X5,m{qz,^) ■ (19) 



There is symmetry between the matrix elements appearing in Eq. (|19|). For example, the only difference between 
Bmrn in the first row and second column and Dmm in the second row and first column is the susceptibility Xj,m{qz, i^), 
with j = 1 or j = 2. Because of the symmetry under the interchange Ri O i?2, the form factor in these two 
matrix elements is the same. Also, the matrix element Cmm may be obtained from Gmm with the replacement 
X3,m{qz,'^) — > Xi.rniqz,'-^) and vice versa. For these two matrix elements, there is symmetry under the interchange 
Ri -!->■ i?3. Further simplification of these results may be achieved for intrasubband excitation (to = to' = 0) by means 
of the identity-2 

d'q±^-—-^Jo{aq±)Joibq±)Joicq±) = h{qza)I(3{qzc)Kn{qzb) 

qi_ + qz 







for a > 5 + c; 6, c, > . (20) 



In the limit when — ^ oo, i.e., when the tubule on the a;-axis is infinitely far away, the inter-tubule Coulomb 
interaction is negligible, and we obtain the following dispersion formula for two tubules on the y-axis, i.e., 

/l + ^,o Co,o \ 

Det 1 + £^0 = . (21) 

V Go^o l + hfi) 

In the limit Oj^ — > oo, the dispersion equation becomes 

/1 + Ao.o Bo,o \ 

Det i?o,o 1 + £^0,0 = . (22) 

V l + /o,o/ 



Thus, as expected, the dispersion equations (21 1 and (22 1 clearly demonstrate a coupling of the modes in a pair of 



the nanotubes leaving the third isolated from the other two in the bundle. We now turn to a consideration of the 
dispersion equation for plasma excitations for a 2D array of nanotubes embedded in a background dielectric medium. 

III. TWO-DIMENSIONAL ARRAY OF NANOTUBES 

Assuming that there is no tunneling between the tubules, the single-particle eigenfunctions for the 2D periodic 
array are 

^ 2 2 

^o.i{.P,z) = e'*'^" e'('=^"-'^-+'=«"''"'')4',fc,,(p-(n,a,e, + nj,ayey)) , 

*.;(P) = ^e"*^<i>,(p) , (23) 

where j = 1,2, •••,M labels the tubules in the nanotube, v = {fc^,, fcj,, fc^} is a composite index for the electron 
eigenstates, ji{p)e^^^'' is the wave function for an electron in the j-ih tubule, with wave vector kz in the axial 
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direction and angular momentum quantum number I = 0, ±1, ±2, • • •, ^f{p) = S{p — Rj), = f^n-x and ky = j^Uy 

with = 0, ±1,±2, • • • and Uy = 0, ±1,±2, • • • ,±^- Here, iV^; = Lx/a^ and Ny — Ly/ay are the numbers of 

nanotubes in the a; and y directions in the array with periodic boundary conditions. Electron motion in the azimuthal 
direction around the tubule is quantized and characterized by the angular momentum quantum number Z, whereas 
motion in the axial z direction is free. Thus, the electron spectrum in each tubule consists of one-dimensional (ID) 
subbands with I serving as a subband index. The spectrum does not depend on kx and has the form 



Consequently, using these results in conjunction with the methods previously employed for the one-dimensional array 
of nanotubes, we obtain the dispersion formula for plasma excitations in a 2D array of nanotubes to be 



Det 



where Gfj — 



2e2 



Jm' (Rj' \J (q. + G^)' + (qy + Glf^ Jra [r, {q. + G^^f + {qy + Gl)' 



qx + G% + i{qy + Gl) 

.^{q. + G%Y + {qy + Giy 



, 



(25) 



and G}(, — 



The result in Eq. ( 25 ) shows that the symmetry of the lattice is maintained in the plasmon spectrum for both the x 



and y directions. In addition, the plasmon excitations depend on the wave vector components q^ and qy with period 
27r 27r 

— and — . In the limit, ay — > cxd, we obtain the following dispersion equation for a linear array of nanotubes on a 

0,^ dy 

2D plane, i.e.. 



Det 



poo 

J2 / dqy 



-Xj'rn'iqz,^^) 



N. 



^J{q. + G%f+qlR,^l J,„ [^^{q. + G^f + qlR^^ 



{qx + G'i,Y + ql+ql 



qx + G% + iqy 
{q,+Gf,y + ql 



, 



(26) 



which agrees with the result obtained by Gumbs and Aizin.^^ 



IV. NUMERICAL RESULTS AND DISCUSSION 



We now turn to numerical calculations for an array of nanotubes as well as a bundle based on the formalism we 
developed in Secs.[TT]and |III| We simulate a metallic graphene tubule in a medium with background dielectric constant 
£s = 2.4 by choosing m* = 0.25me where mg is the free electron mass and Ep — 0.6 eV. The effective Bohr radius 
is aB = h^es/m*e'^ = 1.26 A. All calculations were carried out at zero temperature. We included the transitions 
m — 0, ±1 only in the calculations for plasma excitations in a 2D array. The number of occupied subbands in each 
tubule is determined by its radius, the electron density through the Fermi energy and electron effective mass. The 
occupied subbands are then included in the sum I = 0, ±1, ±2, • • • , zbZmaxi where /max labels the highest occupied 
subband at T = K. For a singlewall cylindrical nanotube of radius 11 A, there are five subbands occupied by 
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electrons corresponding to Z = 0,±1,±2. Lin and Shun^ used the same values of Ss, m* , R and Ep in calculating 
the plasmon excitation spectrum. It was shown that there are three quasi-acoustic plasmon branches associated with 
intrasubband electron excitations with angular momentum transfer m — 0. The reason for this is that the plasmon 
excitations depend on \l\ in this case. There are five optical plasmon branches associated with intersubband electron 
transitions with angular momentum transfer m — ±1. 

In Fig. [2j we present our results for the dispersion relation of intrasubband (to = 0) plasmon excitations of a 2D 
array of singlewall nanotubes. The radius of each tubule is 11 A and the period of the lattice is 35 A. The excitation 
energy is plotted as a function of Qz with = tt/ux and qy — -k j ay. Only the plasmon branches which are not Landau 
damped and which lie outside the single-particle excitation regions are presented in Fig. [2] For the sake of clarity, 
we have omitted the boundaries of the particle-hole continuum in this figure. In Fig. [3] we plot the intrasubband 
plasmon energies as a function of q^ax/'^'K for fixed qy and g^, using the same lattice parameters as Fig. p| Only 
the plasmon modes which are not Landau damped by the single-particle excitations are shown in Fig. [3] There are 
several branches of plasmon excitations in both Figs. |2] and [3] The periodicity of the lattice is preserved in each of 
the plasmon branches when plotted as functions of the transverse wave vector q^. 

It is a simple matter to show from Eq. (26 1 that for a ID array, when g^; = the elements of the determinantal 
matrix with to = 0, to' = ±1 are zero due to antisymmetry of the qj^-integrand. This means that there is decoupling 
of the intrasubband and intersubband excitations. Similarly, it follows from Eq. (25) that when both g^ = and 
qy = 0, the matrix elements of the determinantal equation are zero for to = and m' = ±1 as well. The main result of 
these plots is to show that the Coulomb interaction between the tubules in the 2D array serves to alter substantially 
the dependence of the plasmon frequency on g^. Furthermore, as the separation between tubules is decreased, the 
plasmon frequency is increased due to the role played by the Coulomb interaction between the tubules. Some of the 
plasmon branches in Fig. [3] show a stronger dependence on q^ than others. This variation is more pronounced for the 
high-frequency plasmon mode compared with some of the lower-lying frequency ones. 

The effect of the Coulomb interaction on the intrasubband plasmon excitation spectrum of a 2D array of nanotubes 
consisting of a pair of coaxial tubules is demonstrated in Fig. |4j The outer radius is i?i = 14.4 A and the inner radius 
is i?2 = 11 A. The period of the lattice was chosen as a^^ = = 35 A in both the x and y directions. Only those 
branches which are not Landau damped by single-particle excitations are shown in this figure. The two highest modes 
originate from the high-energy modes of each tubule. The Coulomb interaction between the coaxial tubules leads to 
optical branches which are not present for a single tubule, all of whose modes are acoustic in nature. 

Plasmon excitations arising from intersubband electron transitions for to = ±1 are shown in Figs. [5] and [6] for the 
ID and 2D periodic arrays of single-wall nanotubes. In both figures, only undamped plasmon modes are shown and 
the particle-hole continuum is omitted. The highest modes have a stronger dependence on the transverse wave vector 
qx than the lower modes. This is more so for the ID array than for the 2D array. The main difference between these 
intersubband modes in a single tubule compared with a linear ID array of tubules and when they are arranged in 
a 2D lattice may be explained in the following way. In a single tubule, plasmon modes with to — -1-1 and —1 are 
degeneratd^. The periodic lattice of tubules in either one direction or in two directions removes the axial symmetry 
of the structure. The Coulomb interaction between the tubules lifts the degeneracy and splits each single tubule 
intersubband plasmon mode into two branches increasing the total number of modes to ten for both the ID and 2D 
lattices. The Coulomb interaction seems to have a larger effect on the lowest modes for the ID lattice. Furthermore, 
this splitting is larger for the high-frequency modes compared with the low-frequency ones. This is demonstrated in 
Figs. [5] and |6] As a matter of fact, the splitting of the low-frequency modes is too small to be resolved on the scale 
used in the figures. 

Figure [7|^a) shows the dispersion relation for intrasubband plasmon excitations obtained by solving for the zeros 
of the determinant matrix in Eq. ( |18[ ). The plasmon energy is plotted as a function of the wave vector q^, for fixed 
ttx = dy = 25.0 A. The radius of each of the three cylinders is 11 A. Each of the three cylinders has three plasmon 
branches when it is not interacting with either one of the two other cylinders. The Coulomb interaction between 
electrons on different tubules splits the degenerate plasmon modes for each tubule. The energy gap between split 
modes is very much determined by the distance between the tubules. When ax or ay increases, the separation between 
the split modes decreases, reducing to zero as the limits ax, ay — >■ 00. Our numerical analysis shows that the amount 
by which the branches is split depends on their frequency. To understand the results in Fig. [Tf^a) better, we plot 
the dispersion relation for a pair of singlewall nanotubes in Fig. [T^b) with a^; = 25 A and all other parameters are 
the same as in Fig. [Tf^a). The two highest modes in Fig. [7|^b) are symmetric (ujs) and antisymmetric {(^a) modes. 
When the third nanotube is introduced in Fig. [7][a), a third high-frequency mode appears (a;i,a;2 and Wa) . Since the 
nanotube with its center located at x = a^ is not separated by the same distance from the naotube with its center at 
the origin and the nanotube with its center at y = ay, the plasmon modes uji, uj2 and W3 are not equally spaced. There 
is a total of nine plasmon modes in Fig. and six plasmon modes in Fig. [7][b), but not all of them could be seen 
in these plots. In Figs, [sja) andjsjb), we fixed the value of the longitudinal wave vector at q^ — 0.2 kp and plotted 
the intrasubband plasmon excitation energies as functions of ax- Figure p[a) was obtained for three nanotubes with 
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flj, = 25 A whereas, Fig. |8|b) was calculated for two nanotubes on the s-axis with one of them fixed at the origin. 
The radius of each of the nanotubes was chosen as 11 A and all other parameters are the same as in Figs. [7][a) and 
[7]jb). The highest plasmon modes are degenerate in the limit a^: ^ oo in Fig. [sj^b) for the pair of nanotubes. The uJi 
and UJ2 modes become degeneate as a^; — )■ cx). The mode splits off from these two modes when the Qx coordinate is 
large due to the finite Coulomb interaction between the nanotubes on the y-axis. As a matter of fact, the uj^ branch 
is flat for large Ux since the Coulomb interaction between nanotubes is not affected in this limit. The low-frequency 
branches shown in Figs, [s]^ a) andjsj^b) clearly have a weak depence on the variation of the Coulomb interaction with 
ttx- When Ux ~ 25 A, the uj2 and UJ3 modes anti-cross, indicating that beyond this distance the Coulomb interaction 
is strongest between the two nanotubes on the y-axis. In general, as shown in Eq. (17), the modes with m — 
and m 7^ are coupled to each other. To illustrate the way in which the Coulomb interaction modifies the plasmon 
modes for the three tubules, we only present results for m = 0. In a strict sense, the plasmon excitations cannot be 
categorized as intrasubband and intersubband plasmons as for a single tubule. 



V. SUMMARY 



In concluding this paper, we would like to point out that we used a simple electron gas model to obtain the low- 
frequency plasmon excitation spectrum of a 2D periodic array of parallel nanotubes and a bundle of three nanotubes. 
The nanotubes are assumed dopecP^l in which the charge carriers are introduced onto the graphene tubules by means 
of intercalation, which can be done in carbon fibers or CecP^. Each nanotube consists of coaxial cylindrical tubules. 
The starting point of these calculations was to model the electronic bandstructure of the tubules by a quasi-free 
electron gas confined to the surface of an infinitely long cylinder of finite radius. There was no tunneling of electrons 
between tubules. The random-phase approximation was employed to calculate the plasmon dispersion equation. 
Plasmon excitation energies were obtained numerically for a single-wall nanotube array as a function of the wave 
vector along the nanotube axes, with transverse wave vector Qx- The Coulomb interaction between nanotubes acts 
to split the single-nanotube plasmon modes. We also obtained the periodic dependence of the plasmon energies on qx 
refiecting the translational symmetry of the lattice. Our calculations should serve as a framework for more elaborate 
computations using the tight-binding method for the bandstructure for the higher energy plasmon excitations on 
nanotubes. In our formalism, we showed that we cannot be labeled by the intrasubband or intersubband angular 
momentum quantum numbers due to the Coulomb interaction between charge carriers on different tubules. 
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FIG. 1: (Color on line) Array of nanotubcs whose axes are parallel and taken to be aligned in the z-direction. In 
this paper, we consider a subset forming a bundle, and an infinite periodic array in both x and y directions. 
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FIG. 2: Undamped intrasubband (m = 0) plasmon excitation energy, in units of the Fermi energy Ep, as a function 
of Qz/kp. Here, kp is the Fermi wave number in the ground (l — 0) subband, obtained. The solutions were obtained 
by solving Eq. (25) for q^c = n/ax and Qy — n/ay. The parameters used in the calculation are Ss — 2.4, m* — 0.25me 
where nie is the bare electron mass and ax ^ ay ^ 35.0 A, R — 11.0 A, Ep = 0.6eV. 
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FIG. 3: Undamped intrasubband {m = 0) plasmon excitation energy as a function of the transverse wave vector 
(in units oi 27r/ax)- The solutions are based on Eq. ( pSj ) for Qy — and Qz = 0- All other parameters employed 

obtaining Fig.l2|are the same. 
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FIG. 4: The plasmon excitation energy, in units of the Fermi energy Ep, as a function of Qz/kp, for a 2D array of 
coaxial tubules of outer and inner radius i?i = 14 A and R2 = 11 A, respectively. The lattice constants in both x 
and y directions was chosen as 35 A. The solutions were obtained from Eq. (25), when both a^, — > 00. The 

parameters used in the calculation are the same as in Fig. [2] 
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FIG. 5: For a ID periodic array, the plasmon excitation energy, in units of the Fermi energy Ep, as a function of q^- 
in units of 2-k ja^- Here = 0, the radius of each tubule is 11 A and the period is ax = 35 A. All other parameters 

are the same as in Fig. [21 
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FIG. 6: For a 2D periodic array, the plasmon excitation energy, in units of the Fermi energy Ep, as a function of q. 
in units of 2n/ax- Here Qy — Qz = 0, all other parameters are the same as in Fig. p| 
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FIG. 7: The intrasubband plasmon excitation energy, in units of the Fermi energy Ep, as a function of qz/kp, for 
the three singlewall nanotubes described in Sec. |llj The radius of each tubule was chosen as 11 A and the centers of 
the tubules are at a^; = 25 A and ay — 25 A. in (a). In (b), there are only two nanotubes with one fixed at the origin 
while the position of the second nanotube on the x-axis is at — 25 A. All other parameters are the same as in 

Fig. [2] 
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FIG. 8: The intrasubband plasmon excitation energy, in units of the Fermi energy Ep, as a function of qz/kp, for 

the three singlewaU nanotubes described in Sec.|llj The radius of each tubule was chosen as 11 A and the 
separation between the center of the tubules on the y-axis is ay — 25 A. in (a) while the the coordinate of the third 
tubule on the x-axis is varied. In (b), there are two nanotubes with one fixed at the origin while the center of the 
second one on the x-axis is varied. All other parameters are the same as in Fig. [2] 



